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(57) A method and system are presented for deter- 
mining the best possible quantitative values for one or 
more operational control-variables associated with the 
functioning and control of a nuclear reactor core to iden- 
tify an optimum fuel-cycle design and control strategy 
for operating a particular nuclear reactor core over one 
or more core refueling cycles. The system utilizes a net- 
work of independent processors for contemporaneously 
conducting multiple simulations of a reactor core oper- 
ating under different conditions and constraints repre- 
sented by virtual cases for selected core control-varia- 
bles. The reactor core simulations provide output data 
that is indicative of selected performance parameters 
which reflect the operational state of the reactor 
throughout the duration of a core cycle. As all core sim- 
ulations are completed, the simulation output data for 
each control-variable case is accumulated, normalized 
and mapped to corresponding second order polynomi- 
als which fit the simulator case data for each control- 
variable. Coefficients that uniquely describe each poly- 
nomial are collected by a host processor in an associ- 
ated memory device as a multidimensional data array 
that serves as a type of virtual "response surface". In 
this manner, the response surface acts as a cyber-work- 
space and repository for resultant output data from 
many control-variable case simulations. The polynomi- 
als are then used to predict dependent variable values 
for the performance parameters over a limited range of 
independent control-variable values. The predicted per- 
formance parameter values from each polynomial pre- 
dictor are compared using an objective function to de- 
termine which particular associated independent con- 
trol-variable(s) is likely to produce the greatest improve- 



ment. Another core simulation using the identified val- 
ues is then performed to provide corroboration of the 
polynomial predictors and to calibrate the polynomial 
coefficient data in the response surface with the simu- 
lation process. 
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Description 

[0001] This invention generally concerns nuclear re- 
actor fuel-cycle design and management. More partic- 
ularly, the present invention is directed toward identifi- 
cation of an optimum reactor core design and control 
strategy considering a multiplicity of operational control- 
variables such as fuel composition, fuel exposure, fuel 
placement, control blade positioning and other critical- 
to-quality factors relevant to operation of a nuclear re- 
actor throughout one or more core refueling cycles. 
[0002] A nuclear reactor core includes many individ- 
ual components that have different characteristics that 
may affect a strategy for efficient operation of the core. 
For example, a nuclear reactor core has many, e.g., sev- 
eral hundred, individual fuel assemblies (bundles) that 
have different characteristics and which must be ar- 
ranged within the reactor core or "loaded" so that the 
interaction between fuel bundles satisfies all regulatory 
and reactor design constraints, including governmental 
and customer specified constraints. Similarly, other con- 
trollable elements and factors that affect the reactivity 
and overall efficiency of a reactor core must also be tak- 
en into consideration if one is to design or develop an 
effective control strategy for optimizing the performance 
of a reactor core at a particular reactor plant. Such "op- 
erational controls" (also referred to herein as "independ- 
ent control-variables") include, for example, various 
physical component configurations and controllable op- 
erating conditions that can be individually adjusted or 
set. Besides fuel bundle "loading", other sources of con- 
trol variables include "core flow" or rate of water flow 
through the core, the "exposure" and the "reactivity" or 
interaction between fuel bundles within the core due to 
differences in bundle enrichment, and the "rod pattern" 
or distribution and axial position of control blades within 
the core. As such, each of these operational controls 
constitutes an independent "control-variable" that has a 
measurable effect on the overall performance of the re- 
actor core. Due to the vast number of possible different 
operational values and combinations of values that 
these independent control-variables can assume, it is a 
formidable challenge and a very time consuming task, 
even using known computer-aided methodologies, to 
attempt to analyze and optimize all the individual influ- 
ences on core reactivity and performance. 
[0003] For example, the number of different fuel bun- 
dle configurations possible in the reactor core can be in 
excess of one hundred factorial. Of the many different 
loading pattern possibilities, only a small percentage of 
these configurations will satisfy all of the requisite de- 
sign constraints for a particular reactor plant. In addition, 
only a small percentage of the configurations that satisfy 
all the applicable design constraints are economically 
feasible. Moreover, in addition to satisfying various de- 
sign constraints, since a fuel bundle loading arrange- 
ment ultimately affects the core cycle energy (i.e., the 
amount of energy that the reactor core generates before 



the core needs to be refueled with new fuel elements), 
a particular loading arrangement needs to be selected 
that optimizes the core cycle energy. 
[0004] In order to furnish and maintain the required 
5 energy output, the reactor core is periodically refueled 
with fresh fuel bundles. The duration between one refu- 
eling and the next is commonly referred to as a "fuel- 
cycle" or "core-cycle" of operation and, depending on 
the particular reactor plant, is on the order of twelve to 
10 twenty-four (typically eighteen) months. At the time of 
refueling, typically one third of the least reactive fuel are 
removed from the reactor and the remaining fuel bun- 
dles are repositioned before fresh fuel bundles are add- 
ed. Generally, to improve core cycle energy higher re- 
's activity bundles should be positioned at interior core lo- 
cations. However, such arrangements are not always 
possible to achieve while still satisfying plant-specific 
design constraints. Since each fuel bundle can be load- 
ed at a variety of different locations relative to other bun- 
20 dies, identifying a core loading arrangement that pro- 
duces optimum performance of the core for each fuel- 
cycle presents a complex and computation-intensive 
optimization problem that can be very time consuming 
to solve. 

25 [0005] During the course of a core-cycle, the excess 
energy capability of the core, defined as the excess re- 
activity or "hot excess", is controlled in several ways. 
One technique employs a burnable reactivity inhibitor, 
e.g., Gadolinia, incorporated into the fresh fuel. The 

so quantity of initial burnable inhibitor is determined by de- 
sign constraints and performance characteristics typi- 
cally set by the utility and by the NRC. The burnable in- 
hibitor controls most, but not all, of the excess reactivity. 
Consequently, "control blades" (also referred to herein 

35 as "control rods") - which inhibit reactivity by absorbing 
nuclear emissions - are also used to control excess re- 
activity. Typically, a reactor core contains many such 
control blades which are fit between selected fuel bun- 
dles and are axially positionable within the core. These 

40 control blades assure safe shut down and provide the 
primary mechanism for controlling the maximum power 
peaking factor. 

[0006] The total number of control blades utilized var- 
ies with core size and geometry, and is typically between 
45 50 and 150. The axial position of the control blades (e. 
g., fully inserted, fully withdrawn, or somewhere in be- 
tween) is based on the need to control the excess reac- 
tivity and to meet other operational constraints, such as 
the maximum core power peaking factor. For each con- 
so trol blade, there may be, for example, 25 or more pos- 
sible axial positions and 25 "exposure" (i.e., duration of 
use) steps. Considering symmetry and other require- 
ments that reduce the number of control blades that are 
available for application at any given time, there are 
55 more than 6 million possible combinations of control 
blade positions for even the simplest case. Of these pos- 
sible configurations, only a small fraction satisfy all ap- 
plicable design and safety constraints, and of these, on- 
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ty a small fraction are economical. Moreover, the axial 
-positioning of control blades-aiso-influences the core cy- 
cle energy that any given fuel loading pattern can 
achieve. Since it is desirable to maximize the core-cycle 
energy in order to minimize nuclear fuel cycle costs, de- 
veloping an optimum control blade positioning strategy 
presents another formidable independent control-varia- 
ble optimization problem that must also be taken into 
consideration when attempting to optimize fuel-cycle 
design and management strategies. 
[0007] Traditionally, reactor fuel-cycle design and 
management, including core loading and control blade 
positioning determinations as well as optimization strat- 
egies concerning other variable operational controls, 
are determined on a M trial-and-error" basis based prima- 
rily on the past experiences of the reactor core design 
engineers. Due to circumstances that require a rapid re- 
sponse to changing plant operating conditions, a core 
design engineer may be faced with the formidable chal- 
lenge of specifying values for over 200 independent 
control-variables within a very short time frame. The im- 
pact, for example, of a particular suggested core loading 
arrangement or a control blade positioning arrangement 
on reactorperformanceoverthe duration of acore-cycle 
is usually determined by individual computer simula- 
tions. If a particular design constraint is not satisfied by 
an identified arrangement, then the arrangement is 
modified and another computer simulation is run. Be- 
cause of the relatively long computer simulation time re- 
quired for assessing the impact of a change in the value 
of even a single given independent control-variable, 
man-weeks of human and computer resources are typ- 
ically required before an appropriate fuel-cycle design 
is identified using this procedure. Moreover, using this 
trial-and-error approach, once a fuel-cycle design ar- 
rangement that satisfies all design and safety con- 
straints has been identified, it may turn out that the iden- 
tified arrangement may not provide the actual maximum 
cycle-energy. Therefore, this trial-and-error process 
must continue until the engineers believe that an opti- 
mum fuel-cycle design for the core has been identified. 
In practice, however, it is very possible that a particular 
core arrangement that is not consistent with the engi- 
neers' past experience may be the actual optimum fuel- 
cycle design for the core. Such an actual optimum core 
arrangement, however, may not necessarily be identi- 
fied through the above described trial and error process. 
[0008] Since operational control strategy problems 
generally are considered unique to each reactor plant, 
no known algorithm has provided a viable solution for 
identifying optimum operational control strategies. In 
addition, expert systems have not been applied on a 
broad basis since a standard set of rules typically are 
not really applicable over the wide range of situations 
characteristic of the many different reactor plants and 
types currently in commercial operation. Few method- 
ologies have been developed which can significantly re- 
duce the time required to identify a fuel bundle loading 



arrangement or identify a control blade positioning ar- 
rangement that optimizes cycle energy and satisfies de- 
sign constraints for a wide range of reactors. At least 
one methodology applicable to a wide range of reactors 
5 for identifying optimum control blade positioning ar- 
rangements was developed and is the subject of com- 
monly assigned U.S. Patent No. 5,790,616 to Jackson, 
issued August 4, 1998. 

[0009] Similar methodologies have been developed 
10 for identifying optimum core fuel bundle loading ar- 
rangements. See, for example, commonly assigned U. 
S. Patents No. 5,923,71 7 to Fawks, Jr., issued July 13, 
1999 and No. 5,790,618 to Fawks, Jr., issued August 4, 
1998. For the most part, the above methodologies em- 
's ploy a single processor or computer system to execute 
a specific program routine that simulates the reactor op- 
erating conditions under a selected component arrange- 
ment/configuration of fuel bundle locations or control 
blade axial positions and then the arrangement is opti- 
cs mized by systematically or stochastically evaluating 
possible alternatives. Subsequent to analyzing each po- 
sition or location, random arrangements are created and 
compared with the then best case arrangement identi- 
fied. Another example is a recent methodology for boil- 
25 ing water reactor (BWR) in -core fuel management opti- 
mization that uses a 3-D core simulation computer pro- 
gram called FORMAOSA-B. (See "FORMOSA-B: A 
BWR Incore Fuel Management Optimization Package" 
by B. R. Moore, P. J. Turinski and A. A. Karve, Nuclear 
30 Technology, 126, 153 (1999)). An enhanced version of 
the FORMAOSA-B code has a limited fuel loading pat- 
tern optimization capability through the use of a stochas- 
tic optimization technique called "simulated annealing". 
(See the paper entitled "Effectiveness of BWR Control 
35 Rod Pattern Sampling Capability in the Incore Fuel Man- 
agement Code FORMOSA-B" by A. A. Karve and P. J. 
Turinski presented at the conference on "Mathematics 
and Computation, Reactor Physics and Environmental 
Analysis in Nuclear Application ," published September 
40 1999, SENDA EDITORIAL.S.A.). 

[0010] Unfortunately, the above-described methodol- 
ogies are only applicable for optimizing a single or a few 
operational control-variables at a time such as, for ex- 
ample, fuel bundle loading pattern or control blade po- 
45 sition. Moreover, the above methodologies fail to ad- 
dress the optimization of other important operational 
control-variables such as fuel bundle enrichment, blade 
sequence interval, core water flow, and other independ- 
ent control-variables that may also be critical to quality 
50 and performance. Consequently, it would be highly de- 
sirable to have an efficient optimization methodology 
and apparatus that is broadly applicable over a wide 
range of reactor plant types and which is capable of 
identifying the best possible fuel-cycle core design and 
55 jn-core fuel management strategy in light of all the in- 
trinsic operational control-variables as well as the many 
different specific constraints and considerations that 
may be critical to quality for the operation of a particular 
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reactor plant. 

[0011] An embodiment of the present invention pro- 
vides a system and method for optimizing multiple op- 
erational control-variables of a nuclear reactor to iden- 
tify an optimum fuel-cycle design and develop an oper- 
ational management strategy. In one aspect, the 
present invention is a networked computer system in- 
cluding one or more computers programmed to execute 
a nuclear reactor simulation program and having at least 
one computer programmed to determine the most ap- 
propriate values for selected control-variables that re- 
sult in the optimal physical configuration for operating 
the reactor core over one or more refueling cycles. In 
another aspect, the present invention is a method for 
efficiently determining optimized values for the opera- 
tional control-variables that effect the performance of a 
nuclear reactor core — an operational control-variable 
being the "controllable" physical aspects of the reactor, 
such as fuel bundle loading, control rod pattern, core 
flow, etc., the characteristics of which define the physical 
configuration and operational constraints of a particular 
reactor core. Rather than relying on random or stochas- 
tic search techniques or rule-based techniques in an at- 
tempt to reduce the size of the "search" space, the op- 
timization method of the present invention performs a 
deterministic and exhaustive search for an optimum so- 
lution. 

[0012] In an example embodiment of the invention, a 
plurality of several thousand performance parameters- 
also referred to herein as "dependent" variables be- 
cause of their dependence upon the setting or values of 
the various operational control-variables-are utilized as 
a measure of determining reactor core performance. 
These "performance parameters" include but are not 
limited to parameters conventionally used to gauge re- 
actor core performance, such as CPR, SDM, MAPLH- 
GR, MFLPD, HotExcess, etc. Many of the performance 
parameters analyzed are both spatially and time de- 
pendent, such as, for example, MAPLHGR, MFLPD, 
and MCPR. Consequently, the present invention must 
be capable of analyzing several thousands of such "de- 
pendent" variables. To accomplish this, a preferred em- 
bodiment of the present invention utilizes a plurality of 
processors to conduct separate computer simulations 
covering the operation of the reactor core throughout 
one or more fuel cycles to determine how changes to 
many different control-variables effect the performance 
of the core as measured by the respective performance 
parameters. Preferably, these reactor core simulations 
are conducted utilizing a computer program capable of 
performing three-dimensional (3-D) physics modeling of 
reactor core operation (i.e., the simulator program 
should be capable of handling three-dimensional control 
variables). 

[0013] The method of the present invention signifi- 
cantly decreases the number of required simulations- 
and hence increases the overall computational efficien- 
cy-by mapping the results of a relatively few number of 



reactor core computer simulations to second-order pol- 

ynomialSrThe polynomials are then subsequently used 

to predict quantitative values for performance parame- 
ters (i.e., the dependent variables) over a selected lim- 

5 rted range (i.e., "breadth") of quantitative values or set- 
tings for one or more selected control-variables (i.e., the 
independent variables) . Co nseq u entty, each polynomial 
"predictor" effectively saves the processing time that 
would be required to actually run computer simulations 

10 for the many discrete quantitative values that a particu- 
lar control-variable might assume within a selected 
range or "breadth" of possible control-variable values. 
Each of these polynomial predictors are defined in terms 
of a unique set of coefficient values that are stored in a 

is multidimensional data array in a host computer memory. 
In this manner, the data array serves as a type of virtual 
"response surface" for cataloging and analyzing the re- 
sults of many different simulation cases, based on a 3-D 
physics modeling of the core, in terms of polynomials 

20 indicative of reactor performance that are represented 
. by the polynomial coefficients. 
[0014] Using the polynomial predictors, quantitative 
values for performance parameters are determined for 
discreet control-variable values at selected predeter- 

25 mined increments over the breadth of possible values 
for each control-variable. Each performance parameter 
value is then compared using a conventional "objective 
function", which sets limiting values for each perform- 
ance parameter, to determine the best set of contro l-var- 

30 jable polynomial predictors for optimizing core perform- 
ance. As an option, a particular optimization "resolution" 
level may be selected wherein the combined effect of a 
change in two or more control-variables is analyzed. In 
that instance, predicted values for the two or more per- 

55 formance parameters are combined to generate a net- 
change or "superposition" value indicative of the core 
simulation that would likely result. A corroborating reac- 
tor core simulation is then conducted using the best pre- 
dicted value (s) obtained for each control-variable to pro- 

40 vide corroboration of the polynomial predictors and to 
calibrate the polynomial coefficient data in the response 
surface with the simulation process. 
[0015] The method of the invention presented herein 
can be practiced using most any type of computer net- 

45 work or interconnected system of processors having 
sufficient processing speed and associated data stor- 
age capacity and is not necessarily intended to be lim- 
ited to any particular type of data processor or network. 
Moreover, the software system of the present invention, 

50 including one or more software modules, may be em- 
bodied on a computer-readable medium for transport 
between and/or installation on one or more processors/ 
computers or networked computer systems. In addition, 
the method and system presented herein are believed 

55 to be applicable toward optimizing the fuel-cycle design 
and operation of many different types of reactor plants 
including both boiling water reactors (BWRs) and pres- 
surized-water reactors (PWRs). 
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[0016] The invention will now be described in greater 
— detailrby-way-of example r with- reference to-the draw- 
ings, in which:- 

FIGURE 1 A is a block diagram illustrating a system 
for the optimization of multiple operational control- 
variables for a nuclear reactor; 

FIGURE 1B is a schematic illustration of an exam- 
ple network arrangement of independent proces- 
sors in which the present invention may be embod- 
ied; 

FIGURE 2 is a data flow diagram illustrating the ba- 
sic data flow between processes in an example em- 
bodiment of a software system for implementing the 
reactor core multiple control-variable optimization 
method of the present invention; 

FIGURE 3 is a block diagram illustrating an exam- 
ple embodiment of a software system for imple- 
menting the reactor core multiple control-variable 
optimization method of the present invention; 

FIGURE 4 is a flowchart illustrating functional pro- 
gram control steps performed by a response sur- 
face initialization module in an example software 
embodiment of the present invention; 

FIGURE 5A is a flowchart illustrating functional pro- 
gram control steps performed by a fuel bundle load- 
ing module in an example software embodiment of 
the present invention; 

FIGURE 5B is a flowchart illustrating functional pro- 
gram control steps performed by a control rod axial 
positioning module in an example software embod- 
iment of the present invention; 

FIGU RE 5C is a flowchart illustrating functional pro- 
gram control steps performed by a core flow module 
in an example software embodiment of the present 
invention; 

FIGURE 5D is a flowchart illustrating functional pro- 
gram control steps performed by a sequence inter- 
val module in an example software embodiment of 
the present invention; 

FIGURE 5E is a flowchart illustrating functional pro- 
gram control steps performed by a fuel bundle char- 
acteristics module in an example software embod- 
iment of the present invention; 

FIGURE 6 is a flowchart illustrating functional pro- 
gram control steps performed by a polynomial co- 
efficient development module in an example soft- 
ware embodiment of the present invention; 



FIGURE 7 is a flowchart illustrating functional pro- 
gram control steps performed by a-polynomial us- 
age module in an example software embodiment of 
the present invention; and 

5 

FIGURE 8 is a flowchart illustrating functional pro- 
gram control steps for saving and modifying re- 
sponse surface results in an example software em- 
bodiment of the present invention. 

10 

[0017] The following description is directed toward a 
presently preferred embodiment of the present inven- 
tion, which may be operative as an end-user application 
running, for example, under the Microsoft® Windows 

is 95/NT environment. The present invention, however, is 
not limited to any particular computer system or any par- 
ticular environment. Instead, those skilled in the art will 
find that the system and methods of the present inven- 
tion may be advantageously applied to environments re- 

20 quiring management and/or optimization of any multiple 
control-variable critical industrial/scientific process or 
system, including chemical and mechanical process 
simulation systems, pressurized water reactor simula- 
tion systems, boiling water reactor simulation systems, 

25 and the like. Moreover, the invention may be embodied 
on a variety of different platforms, including UNIX, 
LINUX, Macintosh, Next Step, Open VMS, and the like. 
Therefore, the description of the exemplary embodi- 
ments which follows is for purposes of illustration and 

30 not limitation. 

[0018] Referring first to FIGURE 1 A, a block diagram 
illustrates an example system embodiment of the 
present invention for optimization of multiple operational 
control-variables for a nuclear reactor. Reactor plant 

35 specific design constraints and cycle specific initial data, 
1 , defining a particular reactor core, 3, are provided as 
input data to the optimization system 2 of the present 
invention. Optimized values for operational control var- 
iables (e.g., rod pattern, fuel loading, core flow, etc.) are 

40 provided as outputs for use in the design and manage- 
ment of the nuclear reactor core. 
[0019] Referring to FIGURE 1 B, an example compu- 
ter network arrangement is shown on which the optimi- 
zation method of the present invention may be embod- 

45 jed. A plurality of general purpose computers/proces- 
sors 1 0 are coupled to a local area communications net- 
work (LAN) 15 which may itself be coupled to one or 
more distinct open or private access network(s) 20 for 
communications with one or more remote computers 

50 21 . In a preferred embodiment, the multiple control-var- 
iable optimization method of the present invention is im- 
plemented via software modules resident on at least one 
of computers 1 0. As explained below, the modules may 
be distributed among computers 1 0 or may be resident 

55 on one or more of computers 1 0 (and 21 ) that commu- 
nicate via LAN 15 and/or network(s) 20. 
[0020] As represented in FIGURE 1B, communica- 
tions network 15 and/or 20 can be an open network, 
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such as the Internet, or a private access network, such 
-as a local area network (LAN)-or-a wide-area network 
(WAN). General purpose computers 10 are coupled di- 
rectly or via a modem to network 15 and consist of in- 
dependent processor 1 1 with or without dedicated mem- 
ory 12 in addition to conventional I/O and user interface 
components (not shown). Computers 10 can be any of 
a variety of high speed processors, for example, a VMS- 
Alpha computer system, a Legacy computer system, a 
high-speed work station or a high-speed IBM-compati- 
ble personal computer (such as a desk-top or laptop 
Pentium III processor system). Communications over 
the networks 1 5 and 20 can be accomplished using any 
preferred combination of conventional and proprietary 
protocols that facilitates efficient inter-processor com- 
munications such as, for example, the PCP/IP protocol. 
[0021] Two or more of computers 10 (21), preferably 
systems that are capable of supporting the execution of 
appropriate software for the simulation of nuclear reac- 
tor core operations, are coupled via some communica- 
tions link(s) such as LAN 15 and/or network 20 for ex- 
changing data files and control information. Most any 
conventional reactor core simulation program (or suite 
of programs), such as for example, GE's "PANACEA" 
3-D reactor core simulation program, may be used in 
conjunction with the present invention. This type of sim- 
ulator program is capable of processing three-dimen- 
sional variables defining the core. An input file contain- 
ing values for selected "independent" reactor control- 
variables (e.g., fuel loading, rod pattern, core flow, etc.) 
is provided as an input and the simulator program pro- 
vides an output file comprising values for selected per- 
formance parameters (i.e., reactor core "dependent" 
variables, such as CPR, SDM, MAPLHGR, MFLPD, 
etc.) that are indicative of the state of the reactor core 
at a plurality of discrete intervals (i.e., each and every 
"exposure step") throughout one or more core refueling 
cycles. 

[0022] Referring nowto FIGURE 2, the basic function- 
al processes and data flow within an example software 
system 200 for implementing the multiple control-varia- 
ble optimization method of the present invention are de- 
scribed. Information concerning a selectable "resolu- 
tion" level (explained in greater detail below), other 
processing options and the reactor core cycle-specific 
input data information is preferably input by the user at 
an initial stage (not shown). A cycle-specific reactor core 
profile input file 201 , containing reactor core character- 
istics and operational critical-to-quality constraints spe- 
cific to a particular reactor plant for a particular fuel-cy- 
cle, is built from this user-input information. The cycle- 
specific input data is used to identify initial independent 
control-variable values which define an initial "center- 
point" data case for a particular reactor. This center- 
point data is provided as an input data file 202 to a re- 
actor core simulation program (actual simulation pro- 
gram not shown). A reactor core operation simulation 
207 is conducted using the center-point data. For exam- 



ple, a three-dimensional (3-D) analysis core simulation 

is performed-on-a selected~host"-computer-1 0. When 

the simulation process is complete, a center-point case 
simulation output data file 21 2 is produced. The center- 
5 point case simulation output data from this file is then 
stored in a multidimensional array within the digital stor- 
age memory of the selected "host" computer 10 and is 
used as the basis for creating a type of virtual "response 
surface" 21 9 for evaluating the reactor performance for 
io different control-variable values. 

[0023] Next, separate simulations of the same reactor 
core operating under different physical conditions and 
constraints represented by predetermined changes in 
independent variable values for selected operational 
15 control-variables are conducted contemporaneously by 
the software system. Different simulator input data files 
203-206 are created, each reflecting a change in value 
for a selected control-variable, and each input file is sub- 
mitted to an independent reactor core simulator pro- 
20 gram or process 208-21 1 resident on one or more inde- 
pendent computers or processors (1 0,21 ) connected via 
the communications network (15,20). After performing 
a core simulation based on the values in the received 
input file, each simulator process returns an output data 
25 file 21 3-21 6 reflecting the resultant output values of the 
dependent variables (e.g., CPR, MFLPD, MAPRAT, Hot 
Excess, SDM) of the reactor core. Once all of the reactor 
core simulations for each of the independent variable 
cases 208-211 are complete, the data from simulator 
30 output files 213-21 6 is normalized as indicated at block 
217, for example, by dividing each data item by output 
data obtained from the original "center-point" case 212. 
[0024] After all the simulation case output data is nor- 
malized, the normalized data for each independent con- 
35 trol-variable case is mapped to a set of corresponding 
second-order polynomials. In other words, second-or- 
der polynomials, each of which is characterized by a set 
of associated polynomial coefficients, are selected to fit 
the simulation output data obtained in a few limited 
40 number of reactor core simulations — e.g., three simu- 
lations are preferably used for evaluating each inde- 
pendent control-variable: a center-point case and two 
variation cases wherein the center-point case quantita- 
tive value for the particular control-variable is respec- 
ts tively incremented and decremented. The polynomials 
are then utilized as "predictors" to project quantitative 
values of selected performance parameters for each 
controlvariable. Coefficients which uniquely define each 
polynomial are developed from the normalized simula- 
te tor output data, as indicated at block 21 8, using conven- 
tional algorithmic techniques for solving second-order 
polynomials. This normalized coefficient data is stored 
in an area of computer memory defined herein as the 
"response surface", as represented by block 219. Basi- 
55 cally, response surface 21 9 contains the dependent var- 
iable (performance parameter) "response" of the reactor 
to individual or combined changes in values of the con- 
trol-variables. In this manner, the response surface 
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serves as sort of a cyber- workspace and data-array re- 
— _.pository for-storing4he-resultant-reaGtor core simulation 
output data from different case simulations for multiple 
independent control-variables. 
[0025] Next, all polynomials for each control-variable 
are evaluated 220 and a best polynomial predictor is se- 
lected. As discussed in further detail with respect to the 
Polynomial Optimization And Evaluation Module and 
FIGURE 7, another simulation process 221 is conduct- 
ed using control-variable values provided by the select- 
ed best polynomial predictor to evaluate the modified 
values. If an improvement in reactor performance is in- 
dicated by the simulation results, the modified control- 
variables are accepted as an improvement over the in- 
itial center-point case. This new combination of inde- 
pendent variables is then re-defined as the new center- 
point case and the entire control-variable evaluation 
process is again repeated (as indicated by the dotted 
line in FIGURE 2) until no further significant improve- 
ments are realized. Once it is determined that no further 
improvements are obtainable, the response surface is 
refined using a smaller (more limited) range for control- 
variable values and the above steps are repeated. The 
optimization process as a whole is considered essen- 
tially completed when no further improvements to the 
control-variables are discernible and no feasible reduc- 
tion to the range of control-variable values can be made. 
[0026] In FIGURE 3, an overview of an example soft- 
ware system 300 for implementing the multiple control- 
variable optimization method of the present invention is 
illustrated in terms of functionally related sections or 
"modules" with references to separate accompanying 
FIGURES 4-8 that show example functional program 
control steps for each module in greater detail. One or 
more modules of software system 300, including the 
software system in its entirety, may be embodied on a 
computer-readable medium for ease of distribution and 
installation on one or more processors or networked 
computer systems. Although sections of functionally re- 
lated software are described herein in terms of compo- 
nent software modules that may be individually or col- 
lectively executed by separate processors, the software 
system of the present invention need not necessarily be 
limited to a modular component implementation. As in- 
dicated in FIGURE 3, an example embodiment of soft- 
ware system 300 includes a Response Surface Initiali- 
zation Module 301 , one or more Control-variable Mod- 
ules 302, a Polynomial Coefficient Development Module 
303, a Polynomial Usage Module 304 and a Response 
Surface Save/modify Module 305. A modular arrange- 
ment of the functionally related software within software 
system 300 enhances the overall flexibility and applica- 
bility of the software system to different environments 
by facilitating the use or omission of different Control- 
variable Modules (FIGS. 5A-5E) as desired or appropri- 
ate for a particular application and, moreover, facilitates 
the adding of new and different or updated Control-var- 
iable Modules. 



[0027] Response surface initialization module 301 is 
basically responsible for- accepting-operator-inputted 
data describing operating conditions and constraints for 
a given reactor core (e.g., initial core loading, rod pat- 

s tern, etc.) and creating a starting point or "center-point" 
simulation case for normalizing response surface 219. 
Control-variable modules 302 each contain program 
control steps for creating simulation case data for spe- 
cific types of reactor core control-variables such as, for 

10 example, fuel bundle loading, control rod position, core 
flow, sequence change locations, bundle characteris- 
tics, etc. For each independent control-variable type, 
there may be many independent variable cases to con- 
sider. Furthermore, for each independent variable case 

is considered by a particular control -variable module there 
are at least two core simulations run from which re- 
sponse data is obtained: one simulation is performed 
using the center-point simulation case values with the 
independent variable value increased by a predeter- 

20 mined amount and another simulation is performed us- 
ing the center-point simulation case values with the in- 
dependent variable value decreased by a predeter- 
mined amount. The difference between the increased 
and decreased simulation input values for a particular 

25 control-variable is referred to as the range or "breadth" 
of the control-variable and, since all simulation case re- 
sults are stored in the response surface, it is also re- 
ferred to herein as the "breadth" of the response surface 
(with respect to that control-variable). Each simulation 

30 case result includes the values for all of the operational 
performance parameters (dependent variables) mod- 
eled within the core simulation process. Ultimately, the 
response surface contains at least three core simulation 
case results for each independent variable case: the 

35 center-point case response and two variation case re- 
sponses created by the particular control-variable mod- 
ule. 

[0028] Control-variable modules 302 are preferably 
executed sequentially using a single computer/proces- 

40 sor 10 in the LAN. Additional control-variable modules 
(not shown here) crafted toward particular reactor plant- 
specific considerations may also be used. The control- 
variable modules 302 may be executed in any order and 
any single one or several control-variable modules may 

*5 be used (as indicated by the dotted lines in FIGURE 3) 
depending on the various critical-to-quaiity considera- 
tions and degree of improvement to reactor perform- 
ance that may be desired. Simulator input data files con- 
taining control-variable values are created by each con- 

50 trol-variable module and submitted to other computers/ 
processors in the LAN (or remote network 21 ) that have 
resident core simulator programs. Once a simulation 
case is completed by a processor, it creates a simulator 
output data file containing the resultant values and 

55 sends the file to the computer maintaining the response 
surface. Since reactor core simulations are typically 
very time consuming, this distributed processing ar- 
rangement allows many different core simulation cases 
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to proceed more or less contemporaneously, thereby 
— greatly-reducing-the-overall elapsed-time expended on 
core simulations. Alternatively, different control-variable 
modules could also be resident on different independent 
computers connected within a LAN, WAN or via other 
communications links. For example, in such an embod- 
iment, response surface initialization module 301 resid- 
ing on one computer would place a request over the LAN 
for the execution of a particular desired control-variable 
module to another computer on which that module re- 
sides and then would forward the center-point case data 
from the response surface. 

[0029] Polynomial coefficient development module 
303 contains program control code for mapping the core 
simulation results for each independent variable case to 
unique second-order polynomial curves corresponding 
to each performance parameter (i.e., the operational 
"dependent" variables). The coefficient values of each 
polynomial are determined such that each polynomial 
fits the data from the three simulation cases for its cor- 
responding performance parameter. Polynomial usage 
module 304 contains program control code for exploring 
changes to values of each control-variable, as well as 
changes to combinations of control-variables consid- 
ered together, and determining which changes produce 
the greatest impact on core performance. Since running 
a core simulation is time consuming, the polynomials 
are used as fast predictors (relative to the 3-D simulator 
execution) to determine performance parameter values 
over the input breadth of a control-variable in lieu of run- 
ning a core simulation. The control-variable(s) having 
the greatest performance impact are determined by re- 
iteratively comparing predicted performance parameter 
values against a predetermined objective function. Fi- 
nally, a Save/modify module 305 contains program con- 
trol code for saving and documenting the response sur- 
face and outputting quantified optimum control-variable 
operational values or, alternatively, modifying the re- 
sponse surface if it is determined that results can be fur- 
ther improved by reducing the "breadth" of the response 
surface (explained in greater detail below). 
[0030] Referring now to FIGURE 4, a flow chart illus- 
trates example functional steps performed by response 
surface initialization module 301. The first few initial 
steps 401 -404 basically acquire and identify information 
needed to create an initial center-point simulation case. 
At step 401 , cycle specific reactor core operating con- 
dition data including initial values for control-variables 
(i.e., initial control rod pattern, initial core loading ar- 
rangement, etc.) and an initial response surface breadth 
is specified via operator-input. At step 402, specific op- 
erational constraints which form the design basis of a 
particular reactor plant are identified from the acquired 
operator-input information-such design basis and con- 
straint information aids in the evaluation of an "objective 
function", discussed below, that is used to compare the 
relative quality of alternative solutions. In addition, the 
computer operator may select an input option, dis- 



cussed in greater detail below with respect to the Poly- 
— nomial-Optlmization- And Evaluation-Module and FIG- 
URE 7, that permits the effects on reactor performance 
of a change in the operational value of two or more con- 

5 tro I -variables to be considered in combination. 

[0031] At step 403, the particular independent control- 
variables (core loading, rod pattern, core flow, sequence 
exchange, bundle characteristics, etc.) that are to be 
considered during the optimization are identified based 

10 on the acquired operator-input information. At step 404, 
the fuel bundles to be used within the core are identified 
and sorted according to reactivity value. Next, at step 
405, a core simulation input data file for producing a 
center-point simulation case is generated and submitted 

is to a resident (or remote) core simulation program. Once 
the simulation is finished, the results of the simulation 
are returned in a simulation output file. At step 406, a 
multidimensional array is created in memory as a sim- 
ulation "response surface" and data from the simulation 

20 output file is stored there as an initial center-point case. 
[0032] Next, one or more control-variable modules 
302 are executed to develop simulation case data for 
variations in values for specific control-variables. The 
execution of more than one control-variable module is 

25 optional. Additional control-variable specific modules 
(not disclosed herein) may also be included as desired. 
As previously mentioned, the individual control-variable 
modules may be executed sequentially by a single proc- 
essor or run contemporaneously on different computers 

30 within the LAN or WAN. As the execution of each con- 
trol-variable module results in adding more simulation 
case data to the response surface, the accuracy of the 
present method and the potential reactor performance 
optimization achievable by the present invention is cor- 

35 respondingly enhanced. 

[0033] Referring to FIGURE 5A, the functional steps 
performed by an example control-variable module for fu- 
el bundle loading are discussed first. The fuel bundle 
loading module examines changes in reactor perform- 

40 ance parameters caused by changes in the fuel bundle 
position or loading arrangement. Conventionally, most 
reactor cores are octant-symmetric and, consequently, 
only bundle arrangements within one octant of the core 
need be considered. However, octant symmetry is not 

45 a requirement of the process. As indicated at step 501 , 
it is first determined if fuel bundle loading changes are 
allowed given the pre-identified constraints for the par- 
ticular reactor. If bundle loading changes are not al- 
lowed, program control is passed to another module. If 

so bundle loading changes are allowed, all permissible 
bundle locations are systematically considered by re- 
peating steps 503 through 507 for each different loca- 
tion, as indicated by block 502. 
[0034] At step 503, the known reactivity value of the 

55 bundle at the selected location is changed to a prede- 
termined higher value. A new core simulation input file 
is then generated — the input file reflecting the change 
in fuel bundle reactivity value and a shuffling of the re- 
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maining fuel to minimize any reactivity differences rela- 
-tive-to the center-point.-The core simulation input-file is 
then submitted to an available processor/computer 10 
for simulation processing, as indicated at step 504. (Al- 
though core simulation input files reflecting a "rodded 
depletion" are generally intended and preferred in the 
present example embodiment, non-rodded depletion 
type simulator input files could also be used with the 
present method.) Without waiting for the results of the 
submitted core simulation, the bundle reactivity value for 
the same location is changed, at step 505, to a value 
lower than the original reactivity. 
[0035] The combined amount of increase and de- 
crease exacted to the value for a particular control-var- 
iable, as described herein with respect to the various 
control-variable modules, is predetermined according to 
the particular control-variable being considered and de- 
fines the range or "breadth" of values for which the con- 
trol-variable is examined. 

[0036] Next, at step 506, a new core simulation input 
file having the changed reactivity value is again gener- 
ated and submitted to any available processor/compu- 
ter 10 for processing another simulation. In an example 
embodiment of the present invention, once the simula- 
tion cases in steps 504 and 506 are completed, output 
data parameters from each simulation can be normal- 
ized to the center point, fit to polynomials and stored to 
common response surface 21 9, for example, by each 
processor/computer performing the core simulation. If 
changes in reactivity values for fuel bundles at other lo- 
cations have not yet been simulated, without necessar- 
ily waiting for the core simulations of previous steps to 
complete, a new bundle location is selected and steps 
503-506 are again repeated until all allowable bundle 
locations have been considered, as indicated at step 
507. Ultimately, once all the independent control-varia- 
ble cases for fuel bundle reactivity variations have been 
considered, processing may continue under control of 
another module. 

[0037] FIGURE 5B shows program control steps per- 
formed by an example control-variable module for the 
exploring the different axial positions of the control rods. 
In a manner similar to the fuel bundle loading module of 
FIGURE 5A, two simulation cases for each control rod 
are developed and the simulation results are added to 
the common response surface. At step 509, it is first de- 
termined if control rod pattern changes are allowed giv- 
en the pre-identif ied constraints for the reactor. If control 
rod pattern changes are not allowed, program control is 
passed to another module. If control rod changes are 
allowed, a predetermined control rod is selected for 
analysis, as indicated at step 51 0. Next, at step 511 , the 
initial position value of the selected control rod is in- 
creased by a predetermined amount such that the 
amount of the increase does not violate the physical 
boundaries of the core or the specified user limits. A new 
core simulation input file, having only the selected con- 
trol rod position value changed, is then generated and 



submitted to an available processor/computer 10 for 
- simulation processing, as indicated at-step 51 2. 

[0038] At step 513, the control rod position value for 
the same control rod is changed to a value less than the 

s original position as was done in step 511 . Next at step 
51 4, a new core simulation input file having the changed 
position value is again generated and submitted to an 
available processor/computer 10 for processing a sec- 
ond simulation case. As indicated at step 515, if chang- 

io es in position values for other control rods are to be sim- 
ulated, a new control rod is selected and steps 511-51 4 
are again repeated until ail control rods have been con- 
sidered. As with the fuel bundle loading module, each 
step in control rod positioning module may proceed with- 

15 out necessarily waiting for the core simulations of pre- 
vious steps to complete. Finally, once all the independ- 
ent control-variable cases for control rod position varia- 
tions have been considered, processing may continue 
under control of another module. 

20 [0039] FIGURE 5C shows program control steps per- 
formed by an example control-variable module for de- 
veloping the response surface from changes in the core 
flow. In a manner similar to the other independent con- 
trol-variable modules of FIGURES 5A and 5B, two sim- 

25 ulation cases for each core flow control-variable are de- 
veloped and added to the common response surface. 
At step 519, it is first determined if core flow changes 
are allowed given the pre- identified constraints for the 
reactor. If core flow changes are not allowed, program 

30 control is passed to another module. If core flow chang- 
es are allowed, a particular core flow variable is selected 
for analysis, as indicated at step 520. Next, at step 521 , 
the initial center-point case value of the selected core 
flow variable is increased by a predetermined amount. 

35 a new core simulation input file, having only the selected 
core flow variable value changed, is then generated and 
submitted to an available processor/computer 10 for 
simulation processing, as indicated at step 522. 
[0040] At step 523, the core flow value for the same 

40 core flow variable is changed to a value less than the 
original value similar to step 521 . Next at step 524, a 
new core simulation input file having the changes core 
flow value is again generated and submitted to an avail- 
able processor/computer for processing a second sim- 

4£ ulation case. As indicated at step 525, if changes in core 
flow values for control other core flow variables have not 
yet been simulated, the next independent core flow var- 
iable is selected and steps 521 -524 are again repeated 
until all independent core flow variables have been con- 

50 sidered. As with the other control- variable modules dis- 
cussed above, each step in this module may proceed 
without necessarily waiting for the core simulations of 
previous steps to complete. Finally, once all the inde- 
pendent control-variable cases for core flow variables 

55 have been considered, processing may continue under 
control of another module. 

[0041] FIGURE 5D shows program control steps per- 
formed by an example control-variable module for de- 
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veloping the response surface from changes in se- 
quence interval-ln a manner-similar-to the other control- 
variable modules, two simulation cases for each fuel 
bundle are developed and added to the common re- 
sponse surface 219. At step 529, it is first determined if 
sequence interval changes are allowed given the pre- 
identified constraints for the reactor. If changes are not 
allowed, program control is passed to another module. 
If changes are allowed, a particular sequence interval is 
selected for analysis, as indicated at step 530. Next, at 
step 531 , the initial center-point case sequence interval 
value for the selected sequence interval is increased by 
a user specified amount. A new core simulation input 
file, having only the selected sequence interval value 
changed, is then generated and submitted to an availa- 
ble processor/computer 10 for simulation processing, as 
indicated at step 532. 

[0042] At step 533, the sequence interval value for the 
same fuel bundle is changed to a value less than the 
original value similar to 531 . Next at step 534, a new 
core simulation input file having the changed position 
value is again generated and submitted to an available 
processor/computer for processing a second simulation 
case. As indicated at step 535, if changes in values for 
other sequence interval variables have not yet been 
simulated, a new bundle is selected and steps 531 -534 
are again repeated until ail other relevant independent 
sequence interval variables have been considered. As 
with the othercontrol-variable modules, each step in this 
module may proceed without necessarily waiting for the 
core simulations of previous steps to complete. Finally, 
once ail the independent control-variable cases for the 
sequence interval variables have been considered, 
processing may continue under control of another mod- 
ule. 

[0043] Although the modules depicted in FIGURES 
5A through 5D together demonstrate the ability of the 
present optimization method to consider independent 
control-variables that are capable of having values that 
are considered as "continuous' 1 in nature, such as, for 
example, loading parameters, rod pattern parameters, 
flow parameters, and sequence exchange parameters, 
etc., the method of the present invention can also be 
used to consider changes in "discrete" value control- 
variable, such as bundle characteristics. An example 
control-variable (CV) module for considering discrete- 
value type control-variables is provided using the con- 
text of fuel bundle characteristics as illustrated in FIG- 
URE 5E. 

[0044] Referring now to FIGURE 5E f example pro- 
gram control steps for developing reactor simulation re- 
sponse data from changes in bundle characteristics are 
described. Fuel bundle characteristics, in this example, 
can represent any fuel bundle having differences in axial 
enrichment (e.g., axial Gadolinium variations). Like the 
previously discussed modules, core simulator cases are 
generated and executed for each independent control- 
variable. Upon completion of each independent control- 



variable, the dependent variable output information is 
- -normalized to the- relative-center-point. However, in- 
stead of mapping the response to polynomials, the re- 
sponse is mapped to linear functions. Once all control- 

5 variable modules 302 and corresponding simulation 
cases have finished execution and response surface 
219 has been completed, simulation case data stored 
in the response surface is mapped to polynomials. 
[0045] FIGURE 6 shows example functional program 

io control steps for developing polynomial coefficients for 
mapping each simulation case to a polynomial that fits 
the three data values for each independent variable 
case (i.e., the upper, lower and center-point values). At 
functional step 601 , further processing is delayed until 

15 all of the simulation cases are complete and the re- 
sponse surface has been updated. Next, at steps 602 
and 603, the response surface is accessed and all the 
simulation data produced by control-variable modules 
302 is normalized to the center-point case data. Next, 

20 at functional step 604, coefficients are determined for 
defining a unique second-order polynomial that fits the 
three normalized simulation case values for each inde- 
pendent control-variable. However, since the evaluation 
of certain control-variables (for example, fuel bundle ax- 

25 (at characteristics) can often produce results that are 
much greater than the allowed breadth for that control- 
variable, core simulation results for these type of varia- 
bles are stored in the response surface as discrete first 
order evaluations and are not mapped to polynomials. 

30 Finally, at step 605 , the coefficients for each polynomial 
are saved and further processing continues with the pol- 
ynomial optimization and evaluation module. 
[0046] FIGURE 7 shows example functional program 
control steps for polynomial optimization and evaluation 

35 module 304. This module examines reactor perform- 
ance parameter values predicted by each of the second- 
order polynomials associated with each control-variable 
to determine which control- variable and value produces 
the most significant improvement in reactor perform- 

40 ance. At steps 700 and 701 , polynomials developed 
from each of the control-variable simulation cases are 
accessed from the response surface, sub-grouped and 
used to predict quantitative values for performance pa- 
rameters (e.g., CPR, MFLPD, MAPLHGR, SDM, etc.) 

45 over the breadth of allowable values for that control-var- 
iable. In other words, a control-variable is selected and 
the polynomials associated with each of the perform- 
ance parameters as influenced by that control-variable 
are used to predict a set of performance parameter val- 

50 ues indicative of reactor performance for each of a pre- 
determined number of discrete incremental changes in 
the value of the selected control-variable over the 
breadth (i.e., range of predetermined permissible val- 
ues) of the control-variable. This process is repeated for 

55 every independent control-variable. 

[0047] Under a principle generally known in the art as 
"superposition", the net effect of a plurality of changes 
made to different control-variables together in combina- 
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tion can be determined by the summation of the effects 
of the individual control-variable changes made sepa- 
rately. Accordingly, at the initialization and input stage 
(i.e., when cycle specific inputs and design basis con- 
siderations are identified, e.g. , as discussed above with 5 
respect to steps 401 and 402 of the Initialization Module 
in FIGURE 4), a user of the present system may select 
an optimization "resolution" level as input option that 
permits changes to quantitative operational values for 
more than one independent variable to be evaluated in 
combination with each other. Consequently, if this option 
was previously selected, then, at step 700, the individual 
polynomial-predicted effects of every combination of a 
selected plurality of independent control-variables are 
summarily combined to quantitatively determine the net 
effect that a plurality of changes to different control-var- 
iables made together would have on each of the many 
reactor core performance parameters. The higher the 
selected resolution level, the more independent control- 
variables are evaluated together in combination and, 
hence, the greater the probability of detecting a combi- 
nation that will improve reactor performance. For exam- 
ple, at a selected optimization resolution level of "three", 
a change in the quantitative values for three different 
independent control-variables and every combination of 
three control-variables out of the total number of control- 
variables considered would be evaluated. All discrete 
changes among the plurality of control-variables under 
a particular resolution are examined using the associat- 
ed polynomial predictors for each control-variable. 
[0048] Although higher resolution levels may require 
somewhat longer processing times than lower resolu- 
tion levels, the total processing time is significantly less 
than conventional methods because the polynomial pre- 
dictors are used and combined accordingly instead of 
performing actual computer simulations of the reactor 
core for each case. In this manner, the method of the 
present invention, as exemplified by the present embod- 
iment, is essentially exhaustive and is almost guaran- 
teed to identify the global optimum fuel-cycle design. 
While very high resolution levels may not be feasible in 
practice due to the extended processing time required, 
the capacity of the present invention to permit selection 
of a particular resolution level enables the system user 
to selectively quantify a degree of "closeness" to the true 
absolute optimum which is desired to be achieved. 
[0049] Next, at step 702, for each quantitative value 
change made to a individual control-variable or combi- 
nation of control-variables (i.e., the "independent" vari- 
ables), an "objective function" test is used to quantify 
the relative "worth" or "strength" of that change in terms 
of its effect on improving the performance parameters 
(i.e., the "dependent" variables). The objective function 
sets a particular limiting value for each performance pa- 
rameter that is determined primarily through an integra- 
tion of performance "violations" relative to defined de- 
sign limits, offset by the integration of any performance 
"credits" associated with beneficial results such as ad- 



ditional energy, increased thermal margin, etc. Pre-de- 
termined multipliers (i.e., mathematical factors) are ap- 
plied to design limit values for each of the performance 
parameters-such as, for example, SDM, Hot Excess, 
MFLPD, MAPLHGR, etc.-to provide normalization and 
relative ranking of each parameter. Basically, in step 
702, each predicted performance parameter value is 
tested using an objective function, fashioned in accord- 
ance with conventional knowledge and practice in the 
art, to determine the best set of control-variable polyno- 
mial predictors for optimizing core performance. 
[0050] At step 703, the best values for the control-var- 
iable^) are identified. Since each polynomial predictor 
corresponds to a specific control-variable, polynomial 
predictors are compared, as rated by the objective func- 
tion of step 702, and reiteration of steps 700-702 con- 
tinues until the best values for the control-variable(s) 
have been identified. Next, at step 704, the control-var- 
iable values are compared with the values obtained from 
previous iterations (if any) to determine if any improve- 
ment is found to exist. If no improvement is detected, 
processing continues with the steps shown in FIGURE 
8. If some improvement is found to exist, a core simula- 
tor input case is prepared using the improved values 
from the selected best polynomial predictor(s) corre- 
sponding to one or more control-variables and a core 
simulation is executed, as indicated at step 705. 
[0051] Although the use of polynomials allows for a 
rapid prediction of what changes may constitute an im- 
provement in reactor performance, the core simulation 
at step 705 provides calibration between the simulation 
process and the polynomial coefficient data in the re- 
sponse surface. Essentially, it allows for verifying the 
predicted improvement by providing "actual" (as op- 
posed to "predicted") core simulation data documenting 
the operation of the core under the improved control- 
variables. At step 706, the core simulation results of step 
705 are compared with the core simulation results from 
the center-point case (or the results of previous optimi- 
zations) to determine if any improvement to core per- 
formance has resulted; If the results from the step 705 
core simulation show an improvement over the center- 
point case, the improvement is incorporated and the 
process is repeated again, as indicated at step 708. If 
the results of the core simulation at step 705 is not im- 
proved, the corresponding control-variable(s) is consid- 
ered as "unreliable" and marked as such at step 707. 
Once a predetermined number of unreliable control-var- 
iables is exceeded, as tested at step 709, polynomial 
optimization and evaluation ceases and processing 
continues with the steps shown in FIGURE 8. 
[0052] FIGURE 8 shows example functional program 
control steps for Save/modify Response Surface Mod- 
ule 305. First, the current "breadth" of the response sur- 
face is examined at step 801 (i.e., the breadth of the 
response surface in terms of the range of control- varia- 
ble values explored). If a reduction in the predetermined 
range of values used by the CV modules in creating sim- 
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ulation cases for the control-variables is feasible, then 
-that-range-is-decreased-and-the creation of a new re- 
sponse surface is initiated using the original center-point 
case data. This is indicated at functional step 802 as 
reducing the response surface breadth. At this point, the 
optimization process of the present invention then starts 
over again creating this "new" response surface using 
one or more of the various control-variable modules, as 
indicated by entry point n B" in FIGURE 4. If reducing the 
"breadth" of control-variable values used by the CV 
modules is not feasible, the current response surface 
data is documented (saved) and the optimized control- 
variable values are output, as indicated by steps 803 
and 804. 



3. The method of claim 1 or 2 wherein an objective 
function is-utilized-in step (d) for determining which 

control-variable and quantitative values for a con- 
trol-variable result in optimized reactor perform- 
s ance. 

4. A computer system for optimization of multiple op- 
erational control-variables for a nuclear reactor, 
comprising: 

10 

means for conducting multiple contemporane- 
ous computer simulations of a reactor core op- 
erating under various conditions identified by 
selected quantitative values for a plurality of re- 
's actor core control-variables; 



Claims 

1 . A method of determining an optimum fuel-cycle de- 
sign for a nuclear reactor in terms of quantitative 
values for reactor core control-variables, compris- 
ing the steps of: 

a) conducting multiple contemporaneous com- 
puter simulations of a reactor core operating 
under various conditions identified by selected 
quantitative values for a plurality of reactor core 
control-variables; 



means for mapping selected quantitative val- 
ues of reactor core performance parameters 
obtained from said computer simulations to 
20 second order polynomials; 

means for storing said second order polynomi- 
als in terms of polynomial coefficient values; 

25 means for predicting a plurality of quantitative 

values for performance parameters over a pre- 
determined range of control-variable values us- 
ing said polynomials; and 



b) mapping selected quantitative values of re- so 
actor core performance parameters obtained 
from simulations in step (a) to second order pol- 
ynomials, the selected quantitative values for 
performance parameters comprising values 
obtained from reactor core simulations of step 35 
(a) wherein discrete predetermined changes 

are made to quantitative values of the control- 
variables; 

c) predicting a plurality of quantitative values for 40 
selected performance parameters over a pre- 
determined range of control-variable values us- 
ing the polynomials developed in step (b); and 



means for evaluating said plurality of quantita- 
tive values for performance parameters based 
upon an objective function to determine which 
control-variables and quantitative values for a 
control-variable result in optimized reactor per- 
formance. 

5. A software system having a set of instructions for 
controlling a general purpose computer to deter- 
mine quantitative values for a control-variable opti- 
mized values for multiple operational control-varia- 
bles for a nuclear reactor, comprising: 

a set of instructions formed into a plurality of 
functional modules, comprising: 



d) using an objective function to determine 45 
which control-variables and quantitative values 

for a control-variable result in optimized reactor 
performance based upon the plurality of quan- 
titative values for performance parameters pre- 
dicted in step (c). 50 

2. The method as defined in claim 1 , further compris- 
ing the steps of: 

e) conducting a corroborating computer simu- 55 
lation using quantitative values for control-var- 
iables determined in step (d). 



a) an initialization module having computer pro- 
gram instructions for accepting and storing cy- 
cle-specific input data relating to reactor core 
operating conditions and design constraints; 

c) a plurality of control-variable modules, each 
control-variable module having computer pro- 
gram instructions for initiating computer simu- 
lations of a reactor core operating under spe- 
cific conditions identified by predetermined 
quantitative values of a specific reactor core 
control-variable; 

b) a polynomial coefficient development mod- 
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ule having computer program instructions for 
mapping-second -order -polynomials -from se- 
lected quantitative values of reactor core per- 
formance parameters obtained from reactor 
core simulations, said polynomials being iden- s 
tified in terms of a set of corresponding unique 
polynomial coefficients; and 

d) a polynomial usage module having computer 
program instructions for using an objective 10 
function to evaluate multiple control-variables 
and quantitative values for each control-varia- 
ble to determine optimized values for the mul- 
tiple control-variables, wherein said evaluating 
is based on predicted quantitative values of 
performance parameters obtained from said 
polynomials. 

A software system as set forth by claim 5, further 
comprising a save/modify module having computer 20 
program instructions for saving and/or modifying 
polynomial coefficient data by altering a predeter- 
mined range of control-variable values from which 
reactor core simulations and polynomials are pro- 
duced. 25 

A system for optimization of multiple operational 
control-variables for a nuclear reactor, said system 
comprising a plurality of computer processors inter- 
connected via one or more communication links to so 
form a network, at least one processor acting as a 
host processor having associated data storage 
memory, said memory having stored therein reactor 
core fuel-cycle specific operating condition data 
and reactor core design constraint data, said host 35 
processor programmed to: 

a) initiate a plurality of contemporaneous com- 
puter simulations of a reactor core operating 
under different conditions characterized by pre- *o 
determined quantitative values of one or more 
specific reactor core control-variables; 

b) map second order polynomials from selected 
quantitative values of reactor core performance 
parameters obtained from said computer sim- 
ulations; 

c) determine a plurality of quantitative values 

for performance parameters over a predeter- 50 
mined range of control-variable values using 
said polynomials; and 

d) evaluate said plurality of quantitative values 

for performance parameters based upon an ob- ss 
jective function to determine which correspond- 
ing control-variables and quantitative values for 
a control-variable result in optimized reactor 



performance. 



8. A system for optimization of multiple operational 
control-variables for a nuclear reactor as set forth 
in claim 7, wherein said host processor is further 
programmed to conduct a corroborating reactor 
core computer simulation using optimized quantita- 
tive values for control-variables determined via said 
objective function. 

9. A digital information storage device readable by a 
computer, tangibly embodying a program of instruc- 
tions executable by the computer to perform meth- 
od steps for optimizing a plurality of control-varia- 
bles used in operating a nuclear reactor, said meth- 
od steps comprising: 

a) conducting computer simulations of a reactor 
core operating under various conditions identi- 
fied by selected quantitative values for a plural- 
ity of reactor core control-variables; 

b) mapping selected quantitative values of re- 
actor core performance parameters obtained 
from said computer simulations to second order 
polynomials; 

c) determining a plurality of quantitative values 
for selected performance parameters over a 
predetermined range of control-variable values 
from the polynomials; and 

d) using the plurality of quantitative values for 
selected performance parameters and an ob- 
jective function to evaluate control -variables 
and a plurality of discrete values of control-var- 
iables over said range of control -variable val- 
ues to determine optimized quantitative values 
for said control -variables. 

10. A computer program, embodied on a computer- 
readable medium, for optimizing a plurality of oper- 
ational control-variables for use by a nuclear reac- 
tor, comprising: 

a) computer program code instructions for 
mapping selected quantitative values of reactor 
core performance parameters obtained from 
reactor core computer simulations to second 
order polynomials; 

b) computer program code instructions for pre- 
dicting a plurality of quantitative values for per- 
formance parameters over a predetermined 
range of control -variable values using said pol- 
ynomials; and 

c) computer program code instructions for eval- 
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uating said plurality of quantitative values for 
performance parameters based upon an objec- 
tive function to determine which control-varia- 
bles and quantitative values for a control-vari- 
able result in optimized reactor performance. s 



10 



15 



20 



25 



30 



40 



45 



50 



55 



14 



EP1 113 457 A1 





r 


REACTOR PLANT 




FUEL-CYCLE 


SPECIFIC 




DESIGN & INCORE 


CONSTRAINTS 




FUEL MANAGEMENT 


& CYCLE SPECIFIC 




OPTIMIZATION 


INITIAL DATA 




SYSTEM 



OPTIMIZED 

OPERATIONAL CONTROL 
VARIABLES (ROD PATTERN, 
• • * FUEL LOADING. CORE FLOW, ETC.) 



NUCLEAR REACTOR 
CORE 



Fig.lA 



15 



EP1 113 457 A1 




16 



EP1 113 457 A1 



Fig.2 



REACTOR CORE 
CHARACTERISTIC 
INPUT DATA 
(CYCLE SPECIFIC INPUTS) 



Jl 



■202 



201 



200 



SIMULATOR INPUT 

DATA FOR 
CENTER-POINT CASfl 



■203 



SIMULATOR 
INPUT DATA 



£2 



■207 



Jl 



■204 



-205 



SIMULATOR 
INPUT DATA 



•208 



SIMULATOR 
INPUT DATA 



Jl 



•209 



206 

r 



SIMULATOR 
INPUT DATA 



Jl 



210 



H 



211 



SIMULATOR 
PROCESS 

[XECUTION FOI 

CENTER 
POINT CASE 





SIMULATOR 




SIMULATOR 




r SIMULATOR ' 




r SIMULATOR 




PROCESS 




PROCESS 




PROCESS 




PROCESS 




EXECUTION 




EXECUTION 




EXECUTION 




EXECUTION 




FOR 




FOR 




FOR 




FOR 




INDEP.VARfl 




INDEP.VARfl 


m a • 


INDEP.VARfl 




INDEP.VARfl 




WITH 




WTH 




WITH 




WITH 




INCREASED 




DECREASED 




INCREASED 




DECREASED 




VALUE 




VALUE 




VALUE 




VALUE 


J V A 








^ — - — ' 





212 



SIMULATOR 
OUTPUT 
DATA 



213 



214 



r 



215 



216 





SIMULATOR 




SIMULATOR 




SIMULATOR 




SIMULATOR 




OUTPUT 




OUTPUT 


• • 


OUTPUT 




OUTPUT 




DATA 




OATA 




OATA 




DATA 



( NORMALIZE OUTPUT OATA Y 
T 



•217 



POLYNOMIAL COEFFICIENTS 
DEVELOPED 

I 



r 218 



RESPONSE SURFACE 



i/219 




220 



POLYNOMIAL 

PREDICTOR 

EVALUATION 

AND 
SELECTION 



SIMULATOR 
INPUT OATA 



221 



SIMULATOR 
OUTPUT DATA 



SIMULATOR 
PROCESS 
EXECUTION 



17 



EP 1 113 457 A1 



300 




RS 

INITIALIZATION 
MODULE 

Fig.4 


y30f 








302 


s 

s 

s 


* / \ v 

✓ / \ % 
' / J \ 
/ I \ 

/ I \ 
' I N 




302 

<s 




\ / 


\ ♦ 




FUEL BUNDLE 
LOADING 
MODULE 


CONTROL ROD 
AXIAL POS. 
MODULE 


core aow 

MODULE 


SEQUENCE 
INTERVAL 
MODULE 


BUNDLE 
CHARACTERISTICS 
MODULE 


Fig.5A 


Fig.5B 


Fig.5C 


Fig.5D 
i — j — 


Fig.SE 



X \ 



/ s 





POLYNOMINAL 
COEFFICIENT 
DEVELOPMENT 
MODULE 

Fig.6 


303 






i 






304^ 


POLYNOMINAL 
USAGE 
MODULE 




SAVE/MODIFY 
RS 


Fig.3 


Fig.7 




Fig.8 



y305 



18 



EP1 113 457 A1 



RS INITIALIZATION 
MODULE 

401 



INPUT CYCLE SPECIFIC 

DATA (OPERATING 
CONDITIONS, STARTING ROD 
PATTERN, CORE LOADING, 
RESPONSE SURFACE BREADTH) 



I 



IDENTIFY 
DESIGN BASIS 

(Le. PLANT SPECIFIC CONSTRAINTS) 



IDENTIFY ALL INDEPENDENT 
VARIABLES (LOADING, 
ROD PATTERNS. CORE FLOW. ETC.) 



404 



DETERMINE AND 
SORT BUNDLES BY 
REACTIVITY VALUE 



405 
Vyf 



CORE 
SIMULATOR 



I 



GENERATE INPUT 
FILE FOR CORE SIMULATION 
AND EXECUTE CORE 
OPERATION SIMULATOR 
PROGRAM TO OETERMINE 
INITIAL CASE (CENTER POINT) 



406. 



STORE CENTERPOINT CASE 
RESULTS IN 
RESPONSE SURFACE 



FROM FIG.8 
® 



Fig.4 



(2) TO CV MODULE(S) 
W 5A-5E 



19 



EP 1 113 457 A1 



l~ FUEL BUNDLE 
LOADING MODULE 



1 



CORE 
SIMULATOR 



T 



<& FROM FIG.4 
OR OTHER RS MODULES 




T0FIG.6 
OR OTHER 
RS MODULES', 



~ GENERATE ' "INPUT"F1L£~F0R 

CORE SIMULATION AND SUBMIT/ 1 J 

EXECUTE CORE OPERATION 
SIMULATOR PROGRAM 



-504- 



REPLACE BUNDLE REACTIVITY 
VALUE WITH 

LOWER REACTIVITY VALUE 



( 1 




CORE 




SIMULATOR 


T" — 



505 



GENERATE INPUT FILE 
FOR CORE SIMULATOR 

FILE AND SUBMIT/EXECUTE 

1 



506 



I 

Fig.5A 



20 



EP 1 113 457 A1 



FROM FI6.4 
OR OTHER 
CV MODULES 



CORE 
SIMULATOR 



1 



CORE 
SIMULATOR 



CONTROL ROD 
AXIAL POSITION MODULE 



1 




TO OTHER 
CV MODULES 
OR FIG.6 



51 







YES 


510 


* 


r 




EXAMINE (NEXT) 




CONTROL ROD 




INCREASE VALUE OF CONTROL 
ROD POSITION VARIABLE 



GENERATE INPUT RLE 
FOR CORE SIMULATION 
ANO SUBMI T/EXECUTE 



512 
U 



n 

513 



DECREASE VALUE OF 
CONTROL ROD POSITION 
VARIABLE 



Y 

514 



GENERATE INPUT RLE 
FOR CORE SIMULATION 
AND SUBMIT/EXECUTE 



L 



J 



Fig.5B 



21 



EP1 113 457 A1 



FROM FIG.4 
OR OTHER 
CV MODULES 



CORE 
SIMULATOR 


-+ 








I . 

1 522 



CORE 
SIMULATOR 



519 

ARE 

aow 

-K CHANGES 
.ALLOWED 



CORE FLOW MODULE 



520 



521 




JO OTHER 
CV MODULES 
OR FIG.6 



EXAMINE (NEXT) 
CORE FLOW 
VARIABLE 



INCREASE VALUE OF 
CORE ROW VARIABLE 



± 



523 



J 



GENERATE INPUT FILE 
FOR CORE SIMULATION 
AND SUBMIT/EXECUTE 

I 



-DECREASE- VALUE OF 
CORE ROW VARIABLE 



T 



524 



J"" 1 



GENERATE INPUT FILE 
FOR CORE SIMULATION 
AND SUBMIT/EXECUTE 



J 



Fig.5C 



22 



EP1 113 457 A1 



FROM FIG.4 
OR OTHER 
CV MODULES 



CORE 
SIMULATOR 



CORE 
SIMULATOR 



^7 



SEQUENCE INTERVAL MODULE 

TO OTHER 
CV MODULES 
OR FIG.6 



530 




EXAMINE (NEXT) 
ROD PATTERN 
EXPOSURE 



531 



INCREASE VALUE OF SEQUENCE 
INTERVAL VARIABLE 



532 S 



GENERATE INPUT RLE 
FOR CORE SIMULATION 
AND-SUBMIT/EXECUTL 



J" 



533 



DECREASE VALUE OF 
SEQUENCE INTERVAL VARIABLE 



GENERATE INPUT RLE 
FOR CORE SIMULATION 
AND SUBMIT/EXECUTE 



534 



Z7 



Fig.5D 



23 



EP1 113 457 A1 



r 



BUNDLE CHARACTERISTICS MODULE 



0/? 0r//£>? 

CV MODULES 




EXAMINE (NEXT) 
FUEL 
BUNDLE 



TO OTHER 
CV MODULES 
OR FI6.6 



CORE 
SIMULATOR 



532 J" 



GENERATE INPUT RLE 
FOR CORE SIMULATION 
AND SUBMIT/EXECUTE 



l_ 



J 



Fig.5E 



24 



EP 1 113 457 A1 



r 



POLYNOMIAL COEFFICIENT 
DEVELOPMENT 



FROM LAST 
CV MODULE 

© 



601 



X 



602 



r 



WAIT UNTIL ALL CORE 
SIMULATION CASES 
ARE COMPLETE 

i 



ACCESS SIMULATOR OUTPUT DATA 
FOR ALL RESULTING VALUES OF 
SELECTED DEPENDENT VARIABLES 
(eg.. CPR, MFLPD.MAPLHGR. 
SOM, HOT EXCESS, etc.) FOR 

EACH ROD PATTERN 
EXPOSURE OF EACH INDEP. 
VARIABLE 



603 



NORMALIZE ALL VALUES 
TO CENTERPOINT 



T 



MAP THE NORMALIZED VALUES 
OF--THE-DEPENDENT-VARIABLES— 
j TO A 2nd OROER POLYNOMIAL 
604 S FOR EACH OF THE INDEP. 

VARIABLES (i.e.. SOLVE FOR THE 
COEFFICIENT OF THE POLYNOMIALS) 

; * 

/ SAVE POLYNOMIAL 
/ COEFFICIENTS CORRESPONDING 
. A TO EACH PERFORMANCE 
605 \ PARAMETER FOR EACH 
\ INDEPENDENT VARIABLE 



I 



TO FIG. 7 



Fig.6 



25 



EP1 113 457 A1 



POLYNOMIAL 
OPTIMIZATION I 
AND EVALUATION | 
MODULE, 




26 



EP1 113 457 A1 



©— 

FROM FIG.7 



NO 




SAVE RESPONSE 
SURFACE 



804 



I 



DOCUMENT/OUTPUT 
OPTIMIZED 
VALUES FOR 
INDEPENDENT 
DEPENDENT 
VARIA8LES 



CAN 
'RESPCNST 

SURFACE 
BREADTH BE 
REDUCED, 



802 



,1 



REDUCE RESPONSE 
SURFACE BREADTH 



® 
TO FIG.4 



Fig.8 



27 



EP 1 113 457 A1 



European Patent 
Office 



EUROPEAN SEARCH REPORT 



Application Number 

EP 00 31 1683 



DOCUMENTS CONSIDERED TO BE RELEVANT 



Category 



Citation of document with indication, where appropriate, 
of relevant passages 



Relevant 
to claim 



CLASSIFICATION OF THE 
APPLICATION (tntCI.7) 



D,A 



D,A 



D,A 



DE 198 45 988 A (FRAHATOME SA) 
8 April 1999 (1999-04-08) 

* the whole document * 

DE 44 16 463 A (ASEA ATOM AB) 
17 November 1994 (1994-11-17) 

* the whole document * 

PATENT ABSTRACTS OF JAPAN 

vol. 016, no. 402 (P-1409), 

25 August 1992 (1992-08-25) 

& JP 04 132995 A (POWER REACTOR & NUCLEAR 

FUEL DEV CORP), 7 May 1992 (1992-05-07) 

* abstract * 

EP 0 786 782 A (GEN ELECTRIC) 
30 July 1997 (1997-07-30) 

* the whole document * 

EP 0 823 712 A (GEN ELECTRIC) 
11 February 1998 (1998-02-11) 

* the whole document * 

US 5 790 618 A (FAWKS JR JAMES E) 
4 August 1998 (1998-08-04) 

* the whole document * 



1-10 



1-10 



1-10 



G21D3/08 



1-10 



1-10 



1-10 



TECHNICAL FIELDS 
SEARCHED (lnt.CI.7) 



G21D 



The present search report has been drawn up tor afl claims 



THE HAGUE 



Del* ol completion & tt* Mart 

17 April 2001 



Deroubalx, P 



CATEGORY OF CITED DOCUMENTS 

X : pa rt feu tarty retomuit it taken alone 

Y : particularly relevant if combined with another 

document of (he same category 
A : technotogcai background 
O : non-wrtnen disclosure 
P: 



T : meory or principle underlying the hvemlon 
E : earlier paten) documeni. but pubtshed on. or 

after Ihe filing date 
O : documeni cfted in the application 
I. : document cited lor other reasons 

A : member of the same patent family, corresponding 



28 



EP1 113 457 A1 



ANNEX TO THE EUROPEAN SEARCH REPORT 
ON EUROPEAN PATENT APPLICATION NO. 



EP 00 31 1683 



This annex lists (he patent (amity members relating to the patent documents cited in the above-mentioned European search report. 
The members are as contained in the European Patent Office EDP file on 

The European Patent Office Is in no way liable for these particulars which are merely given for the purpose of information. 

17-04-2001 



Patent document 
cited in search report 



Publication 
date 



Patent family 
member(s) 



Publication 



8 



DE 


19845988 


A 


08-04-1999 


FR 


2769402 A 


09-04-1999 










BE 


1011957 A 


07-03-2000 


DE 


4416463 


A 


17-11-1994 


SE 


509235 C 


21-12-1998 










JP 


6347586 A 


22-12-1994 










SE 


9301613 A 


12-11-1994 


JP 


04132995 


A 


07-05-1992 


NONE 






EP 


0786782 


A 


30-07-1997 


US 


5923717 A 


13-07-1999 










FI 


970353 A 


30-07-1997 










JP 


9304570 A 


28-11-1997 


EP 


0823712 


A 


11-02-1998 


US 


5790616 A 


04-08-1998 










JP 


10123284 A 


15-05-1998 


US 


5790618 


A 


04-08-1998 


NONE 







u For more details about this annex : see Official Journal of the European Patent Office, No. 12/82 



29 



